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ABSTRACT 


Large, fresh-water ice blocks (0.5m X 1.25m X 5m) were 
towed in Monterey Bay at speeds in range from 0.7 to 1.2 knots. 
The objectives of the experiments included measurement of 
gross regression rates of ice surfaces, wake temperature, 
turbulent thermal boundary layer, ice interior temperature 
profile, and observation of shape changes over the melting 
period. The research was conducted over a period from 
October 1977 to January 1979. The measured regression rates 
at several points on the ice blocks were compared to theoret- 
ical predictions uSing a turbulent flow ice ablation model 
developed by Dr. Owen Griffin of the Naval Research Laboratory. 
Griffin's model predicted a regression rate of 280mm/hr 
compared with the measured value of 260mm/hr at a selected 
point. Temperature profiles in the ice interior and ice-water 


boundary layer-were used in a one-dimensional energy balance 


across the melting interface. Ice ripples were observed on 
all ice block tests. The ripples are important for both heat 
transfer and drag considerations. The wavelengths of the 


ripples were measured and the average wavelength determined. 
The measured average wavelength was 6.1 inches compared to 
the theoretical average predicted by Tatinclaux and Kennedy 


of 5.8 inches for a selected experiment. 
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Zi. INTRODUCTION 


eeoUotlTrFICATION FOR STUDY 

There has been Significant interest in recent years 
concerning the movement and deterioration of the large ice 
bodies that drift from the regions of the Arctic and 
Antarctic to the warmer waters oeene Northern and Southern 
Hemisphere. The United States Navy in particular has been 
concerned about the impact these icebergs have on the 
temperature, density, and salinity of surrounding waters. 
Changes in these characteristics result in Significant 
modification of the speed of sound through the water, which 
in turn complicates the problem of submarine detection, 
Submarine navigation, and bottom navigation. The Navy and 
Coast Guard want the capability to predict with some degree 
of certainty the "life" of an iceberg and thus how it will 
affect surface ship operations over a period of time. In 
addition, as the need for additional sources of crude oil 
Seeeands, Offshore drilling in near-arctic regions surely 
will necessitate the movement of an iceberg by artificial 
@ems CO avoid drifting into a drilling rig. How far to 
push or pull the iceberg, the power required, and the 
effect on the iceberg as a result of this movement are 
questions that must be dealt with before such an undertaking. 
A most important interest in iceberg deterioration, is 


the possibility of supplying water from icebergs to the 


Vee 





extremely arid regions of the world. A discussion of this 


area is given in Appendix A of this thesis. 
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MDC WwohOUND OF STUDY 


This thesis concerns a study conducted by the Naval 
Postgraduate School (NPS) in conjunction with the Naval 
Research Laboratory (NRL) on the ablation of ice bodies 
moving thorugh seawater. Relatively large (16' x4' x1.5') 
ice models were manufactured and moved through a body of 
Salt water (Monterey Bay, California). In addition, smaller 
models were placed in a laboratory flow channel for experi- 
ments and analysis. The ultimate goal was the ability to 
predict the overall ablation rate of an ice body moving 
through seawater. The prediction should be verified by 
comparing the experimental results for this thesis with 
ablation theory previously published by Owen Griffin [1,2] 
or NRL. 

A first important step in the study was to analyze the 
scaling problem relating a large iceberg (15,000' x 2500' x 700°) 
to a much smaller model. For the tests in Monterey Bay, 
the size of the model was constrained by manufacturing, 
crane, and transportation facilities available. In addition, 
for the experiments in the flow channel, the model size 
was limited by the flow channel equipment at NPS. The 
scaling 1S important because the experimental results are 
more valuable if the information can be extrapolated to a 
full-scale iceberg. The physical effects of the ice moving 


relative to the seawater (i.e.: shape change, ripple 


ge, 





effect, and wake temperatures) were examined. The thermal 
boundary layer between the ice and water was measured and 
compared with that theoretically predicted. Finally the 
heat transfer through the ice was analyzed by comparing 
temperature profiles in the ice over a period of time with 
those predicted from a solution of the equations governing 


transient heat transfer. 
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mir. ICE MODEL SCALING AND INITIAL LABORATORY TESTS 


A. GENERAL 

The problems encountered in modeling a large ship for 
laboratory water channel experiments are also inherent, 
but to a greater magnitude, in modeling a large ice body. 
This is due to several factors, the most predominant being 
the immense size of the ice bodies to be modeled coupled 
with the low speed of transport (1.0 kt). Also resolving 
the hydrodynamic requirements with those of the heat 
transfer analysis to realize an acceptable model from which 
melt rates can be extrapolated to the full size ice body is 
also an important consideration. Water channel and towing 
tank modeling in general are discussed thoroughly in Chapter 
foe ref {3], and modeling particular to ice bodies is 
Seeeussed by Murphy [4] and Griffin [1,2]. Murphy dis- 
cussed the methods used to determine towing forces and mass 
loss during transit, and Griffin treated the heat transfer 
PLoplem specifically to predict ablation rates of large 


lce bodies. 


Pees DRODYNAMIC REQUIREMENTS 
Two well-understood dimensionless groups of variables 
are commonly used to relate the forces on a body moving 
through a fluid. These are: 
VL inertial effects 


ee > Bumber (we) = ) ~ viscous effects ‘*’ 


Za 





Froude number (Fn) = = = SS SES (2) 
The Reynolds number is important in defining whether a 
meme OVer a body is laminar or turbulent, a fact that has 
implications in both towing force and heat transfer analyses. 
Since the kinematic viscosity (v) is a property of the 
fluid, the only variables for Reynolds number which can be 
easily controlled by the experimenter in ice model testing 
are the velocity of the fluid (v) and the characteristic 
length of the model (L). For ideal modeling, the Reynolds 
numbers of the model and prototype should be identical to 
assure the same laminar or turbulent flow conditions. This 
means that a much smaller model must move through the fluid 
at a very high velocity. The high speed of the test fluid, 
however, iS in direct conflict with the Froude number which 
1s used in determining residuary resistance of a prototype. 
In addition, the Froude numbers of the model and prototype 
should be dagen for ideal modeling. Note that the 
Froude number also has only fluid velocity and model length 
aS independent parameters. 

The solution commonly used in ship modeling for the 
conflicting Reynolds number - Froude number requirements 
is to scale the velocity and model size to satisfy Froude 
number similitude. Reynolds number scaling is unavoidably 
violated. However, models are tested at either moderately 


high Reynolds numbers (10°) to assure natural transition to 


Ze 





turbulent flow over forward regions of the model or with 
the use of artificial boundary layer trip mechanisms placed 
so as to stimulate transition at the predicted natural 


transition point on the model. 


SeeeeHeAL TRANSFER ANALYSIS REQUIREMENTS 

Griffin [1,2] at NRL proposed a heat transfer model 
for a large glacial ice body moving through seawater with 
turbulent flow. In Griffin's model the local turbulent 


melting rate is given by: 





STURB DSty TURB (3) 
where: 
STURB Local turbulent melting rate (m/s) 
D = Stefan number parameter determined by: 
A 
a 1 +B i) 
and: 
~-l 
A = Seawater Stefan number = Cop (Tt - T2)M 
B —eeECe Steran number ~- C__(T.{ - T ym t 
PS ' F Oo 
= (5) 
Sty TURB Local turbulent Stanton number 
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where Ris the ratio of interfacial heat transfer with 
and without melting; R is assumed to have the same value 
in laminar as well as turbulent flow. 

The local turbulent Stanton number on a flat plate 
1s: 


Nu 0297 Re “1/9 
x pps 


St = = Sas ee 
h, REF Re, Pr 1+1.48 Re, 1/10 5-1/6 (p,-4) 





Peyesubstituting eq. (6) into eq. (5), and then into eg. 
(3), the variables affecting the local turbulent heat 


transfer rate can be seen as: 


0297 Re 2/9 
a = i gt (7) 
TURE 1+#1.48 Re “1/10 (Pr-1) 


Since Prandtl number (Pr) and Stefan number (D) are 
properties of the seawater and ice, both numbers are the 
same for model and prototype. Likewise R, which is depen- 
dent on the properties of ice and the solid-liquid phase 
transformation, is also invariant between model and 
PmeEOcype. 


For ideal modeling: 


SY 


eS = 1 (8) 
“prototype 
From eq. (7) it is obvious that agreement between Reynolds 


numbers of the model and prototype is paramount for ideal 


modeling as stated by eq. (8). 


24 





From the above consideration, model testing of ice 
bodies is limited to a regime of Reynolds numbers and 
Froude numbers as shown in Figure l. The points indicated 
on the figure are for experiments that will be discussed 


later in this chapter and in succeeding chapters. 


D. WATER CHANNEL FACILITIES AT THE NAVAL POSTGRADUATE SCHOOL 

The Naval Postgraduate School. has several water channels 
which are pictured in Figures 2-5. The characteristics of 
each channel are summarized in Table I. Since buoyancy was 
considered an important factor in the heat transfer modeling 
analysis, the water tunnel shown in Figure 2 was eliminated 
as a testing facility. The water channel in Figure 3 was 
also eliminated because of the extensive modifications 
required to obtain a uniform velocity profile over the 
width of the channel. The remaining two water channels 
were utilized; however, their relatively small size pro- 
hibited the achievement of Reynolcs numbers and accompanying 
turbulent flow conditions equivalent to those on a large 
ice body. The inability to achieve these conditions was 
a direct result of the conflicting requirements of high 
water channel velocities (2-4 ft per sec) and long model 
lengths (3-6 ft) necessary for moderately high Reynolds 
numbers and the low velocities needed for satisfactory 
Froude number scaling. 

Large water channels, Similar to the new Environmental 


Protection Agency water channel in North Carolina, with 
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Beemensions of L = 80 £t, W= 8 ft, and D = 4 ft, would 
allow the testing. of models large enough to achieve a more 
Satisfactory scaling of Reynolds number, geometric size, and 
velocity. This scaling would be achieved by virtue of the 
relatively long models (8-16 ft) that could be tested at 
the relatively low velocities (0.1-0.5 ft per sec) required. 
However, the time, expense, and coordination involved in 
conducting experiments remote ae NPS were deemed prohibi- 
tive; therefore the decision was made to conduct a series 
of large scale model tests in Monterey Bay adjacent to 
NPS. The details of the experiments in the bay are dis- 
cussed in Chapters IV-VI. 

A review of literature dealing with large ice bodies 
revealed no experimental data on melt rate, shape change, 
Or regression rates of fresh water ice moving through sea- 
Water. Therefore, in order to gain at least a qualitative 
understanding of what these characteristics were, a series 


of model tests were conducted in the NPS water channels. 


fee LNITIAL LABORATORY EXPERIMENTS 

The water channel selected for the initial tests is 
pictured in Figure 4 with characteristics listed in Table I. 
The size of the ice models was selected to avoid both 
blockage and side and bottom wall flow interference in 
the channel. 

Four rectangular pieces of ice were constructed with 


= 


ee eci mensions Of L = 8.0 in, W= 1.5 in, and D = 1.3 in, 
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Figure No. 


TABLE I 


Naval Postgraduate School Water 
Channel Characteristics 


Test Section Velocity 
Dimensions Range 
Eve oH (in) ft/sec 
16/4/8 0 - 30 
96/48/5 0-0.5 
Soy 2ZyrG O32 ==320 
36/24/8 eS 1020 


2a 
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thus achieving a satisfactory ratio of channel width to 
Megdemewidth of 10:1. The models were frozen in a standard 
freezing compartment at a temperature of 10°F. A wire 
restraining device was frozen into the body of the model for 
handling during testing. To test the effect of a stream- 
lined shape during identical flow channel conditions, one 
model was constructed with all sides sloping at 30 degrees 
from the vertical. Overall dimensions at the top surface 
fee, = 9.7 in, W = 3.7 in, and D= 1.2 in. The ratio of 
channel width to model width for this model was a less 
Satisfactory 4:1; however, the low channel water flow 
velocity virtually precluded any significant wall inter- 
ference. 

Each model test was similar and consisted of transferring 
the model from the freezing compartment to a scale immedi- 
ately above the flow channel where the model was weighed 
and photographed simultaneously in the side view. The model 
was lowered inot the flowing water and allowed to float 
naturally. Directional control of the model was accomplished 
by attaching the model restraining device to a cross-beam 
Over the channel. At regular intervals (1 min), the model 
was removed, weighed, and photographed until the restraining 
device separated from the model. Dyes were frozen into the 
models, or, alternatively, were injected in front of the 
meee! fOr visualization of the flow patters. In four out 


of five tests conducted, the data on weights of the receding 
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Mocels, air and water temperature, flow velocity, and 
regression rates at several areas of the models were noted. 
The total time in water for each model was as follows: 
meeezemin/50 sec, (2) 2 min/30 sec, (3) 3 min/30 sec, 


feeeinvalid (due to hookup difficulties), and (5) 2 min/30 sec. 


PeeeewAlhR CHANNEL TEST RESULTS 

The first four models melted through a series of shapes 
Semibbar to those illustrated in eee 6. The second model 
tested 1S shown in Figure 6; results of the melting are 
given in Table II. The rates of regression at three 
Selected surfaces for the model in Figure 6 are given in 
iaote Tit. 

The frames of Figure 6 are side views with the front 
Meee 1Ce Orlented to the right of the photograph. The 
streamlined ice model was the fifth to be tested. The 
sequenced melting shape configuration for this model is 
Shown in Figure 7 with similar orientation. It was noted 
that despite the significant geometric differences between 
the streamlined model and the four previous models, the 
melting configuration is practically identical. The regres- 
Sion on the front, bottom, and rear of the models can be 
seen in Figures 6 and 7 which clearly illustrate the area 
of highest regression as that on the bottom, approximately 
three-tenths of the model length from the front. The 
laboratory conditions included relatively small Reynolds 


numbers (approx. 10°) compared with the large numbers 
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TABLE II 
Melt of-Second Rectangular Ice Model 


(See Figure 6) 


Photo Frame Weight (g) Time (sec) 
Upper ie 2 0 
Upper Middle 2 60 
Lower Middle 49 120 
Lower is: Se 

mee Temperature - 23°C Water Temperature 21.3°C 


Hlow Volocity — 0.95. ££/sec 


TABLE III 
Regression Rates of Second Rectangular Ice Model 
(See Figure 6) 


Average Rates for 120 Seconds 


Rear Bottom (neck) Patea@ mde 
ao OnahiGg. as fetoe 
sec sec sec 
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(approx. 107) associated with large Atlantic ice bodies. 
Given these small Reynolds numbers eee in the labora- 
tory, it was concluded that the flow conditions in the 
laboratory were largely laminar and that the decisive necking 
effect observed was associated with the reattachment of the 
laminar boundary layer initially separated from the model 
by the sharp corners of the blunt bow. 

For the melting observable in Figures 6 and 7 along the 
profile surfaces of the models, the Nusselt numbers, which 
measure the effectiveness of heat transfer, can be calcu- 
lated to show this effectiveness at any distance along the 
profile surface. Figure 8 shows a plot of Nusselt number 
memaerunction of length along the ice model. The plot 
Shows a well-defined region of intense heat transfer at a 
distance of approximately 2D (D = Depth) from the front of 
the model and an accompanying decay of melting at approximately 
3.5D from the same point. A representative theory for 
laminar flow on a flat plate with constant surface tempera- 
rer is also shown on the plot in contrast to highlight the 
empirical results observed in the laboratory. 

To explain the high melting rate in the forward third 
region on the models, the following hypothesis is proposed: 
Due to the initially sharp corner on the model, the flow 
separates from the model as it passes the corner point. 

In the region 0-2D, there occurs a recirculating flow with 


reattachment near X = 2D. As is well known, the area where 
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Beatetachment occurs is a region of high heat flux as analyzed 
by Meyer [5] and Koseman [6]. Figure 9 illustrates this 
hypothesis. 

An alternate hypothesis involves a trapped vortex 
filament as shown in Figure 10. The sense of the vorticity 
is shown by the arrows and follows the customary right-hand 
rule. The vorticity originates because of the flow separa- 
mem from the blunt bow. ini cuey the vorticity is shed; 
but as the ice ablates, the vortex filament becomes trapped. 
Due to the vorticity, the warm channel water 1s circulated 
close to the ice resulting in large temperature gradients 
and accompanying high heat flux. The observations made during 
the model experiments were consistent with either the 


reattachment or trapped vortex filament hypothesis. 


G. CONCLUSIONS 

The normal problem of resolving the conflicting require- 
ments of Reynolds number and Froude number scaling in 
modeling ships is magnified when applying standard modeling 
techniques to large ice bodies. This difficulty in modeling 
is primarily due to the extremely large size and low trans- 
port speed of these full-size bodies. The high Reynolds 
number required for the model is necessary to assure turbu- 
lent flow and to minimize the difference in theoretical 
melt rate between the small-scale ice model and full-size 
Meensbody. The only practical controlling variables available 


in developing the model are the size of the model and the 
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velocity of the water. The flow separation around the 
frontal corners of small models results in an extremely 

high heat flux and accelerated melting rate at the point 

of reattachment. This severe melting tends to mask all 
other characteristics that could be determined from such 

a small model. The test results indicated that large scale 
model testing in Monterey Bay, despite the lesser controlla- 
bility of environmental cones are needed in order to 
determine more accurately the melt-rates, shape changes, 

and other characteristics of a full-scale ice body underway 


in seawater. 
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IV. LARGE-SCALE ICE MODEL EXPERIMENTS 
IN MONTEREY BAY 
The following is a discussion of the ice model size 
and shape determination, freezing container design, freezing 
technique, and general test procedure and a summary of 
each test conducted in the Bay. Specific data results 
and theoretical correlations from the various tests are 


presented in Chapter VI. 


eo 4e AND SHAPE DETERMINATION OF MODEL 

Several major considerations of both scientific and 
practical nature were important in determining the minimum 
and maximum size of the ice body model to be tested in 
Monterey Bay. The controlling minimum dimension was the 
model length, which must be sufficient to achieve natural 
transition to turbulent flow at speeds of 1-2 ft/sec (Re 
Or). The model required sufficient volume to withstand 
the following test conditions: a one-hour loading and 
transit time (20 mi) from the vendor to the launch site on 
an open truck with minimal insulation, a second hour period 
at the launch site for testing preparations and launch, a 
30-minute underway period in the Bay, and recovery from the 
Meter intact for examination of resulting changes in physical 
characteristics. The maximum size of the model was limited 
primarily by the size of the freezing compartment which 


was the only large compartment available in the immediate 
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geographical area. However, other factors such as the 
load capacity of a flatbed semi-truck, eieons Side space 
restrictions, the available crane lifting capability, and 
the maneuverability of a boat with ice in tow within the 
restricted harbor were also important in determining the 
maximum size. The utmost consideration throughout the 
development of this previously untried test series was the 
complete safety of all personnel during the tests. 

The model shape was desired to simulate a tabular 
Antarctic iceberg similar to those that continually break 
off from the ice shelves and drift Northward. These ice- 
bergs have considerable length as compared to width and 
depth and thus have good longitudinal hydrodynamic stability 
for towing. The ends of these icebergs are blunt and have 
a flat table-top appearance above the waterline. The 
tabular iceberas are readily duplicated in shape using an 
easily manufactured rectangular container. 

The result of the foregoing size and shape considera- 
tions was an ice model with dimension of L = 16 ft, W= 4 
Mmemand © = 1.5 ft. This particular size and shape satis- 
fied all of the requirements and was used unchanged through- 
Out the entire series of experiments conducted in Monterey 


may . 


BeeeecONSTRUCTION OF THE ICE MODEL FREEZING CONTAINER 
The container was manufactured by a local vendor from 


galvanized steel sheeting with welded seams. The 
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specifications for the container are listed in Table IV, 

and the container iS pictured in Figures 1l and 12. Suffi- 
cient integral bottom support structure was included to 
facilitate the lifting of the 6000 lb. ice model and test 
accessories. Removable plugs were installed to allow flooding 
of seawater into the container during launching and draining 
of this water after launch of the model. This simple design 
proved successful requiring no modifications throughout the 


test series. 


fencer MODEL CONSTRUCTION 

For the first two experiments conducted in Monterey Bay, 
the container was placed in the vendor's freezing compartment, 
filled to the 18-inch level with demineralized water (18 ppm 
dissolved solids), and frozen as a single layer. Prior to 
fering with water, the towing bollard, seen in Figure 13, 
was suspended over the container centerline, 48 inches 
behind the front wall, with the bollard base 9 inches above 
[m@emooctom for frrezing into the ice. The bollard position 
was determined to afford maximum lateral towing stability, 
while remaining intact in the ice during the course of the 
test with rapid melting of the front and bottom surfaces 
of the ice. This position also precluded any hydrodynamic 
interference with the water flowing past the ice and the 
Subsequent masking of natural physical characteristics by 


mags Llow interference. 
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TABLE IV 


Spcciplcartons OLfeice Model Container 


Dimensions (bottom) 
iLergkencle, iia aes 
ect = 4 Ft 


Heng etl 66. Ft 


Slope of Sides 


Seven degrees out on all sides from bottom. 


Material 
Sides - 12 GA galvanized mild steel. 
Bottom support structure - mild steel of 


sufficient strength to support 6500 lbs. 


Drain Plugs 
imseal! two 2° sdlia drain plugs on Center line 


BOEEOM OL ErOnt and rear sides. 
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Mae Size Of the bollard vertical post and the horizontal 
base plate were designed to minimize the towing stress 
concentrations on the ice. This was accomplished by dis- 
tributing the towing forces over a wide area in both 
Memezontal and vertical planes, without significantly 
weakening the ice at the critical vertical lateral lane 
passing through the bollard installation point. The bollard 
protruded above the ice, as pictured in Figure 12, and also 
served as an anchoring point for instrumentation during the 
fEescs. 

Two problems developed from freezing the ice as a Single 
Mmeemen layer. First, the expansion of the ice during 
freezing resulted in a large upheaval of the top surface 
as the core solidified. The upheaval of the ice was accom- 
panied by large vertical fissures up to one inch in width 
on the surface and extending several inches into the ice 
at some points as shown in Figure 12. Second, as a result 
of the upheaval, the effective freeboard (distance from top 
surface to waterline) was Significantly lowered allowing 
seawater to flow over the top surface of the ice. The result 
was rapid erosion of the upper surface of the bow and 
fissured areas with premature disintegration of the model 
as shown in Figure 14. 

To eliminate the model's top surface upheaval and low 
freeboard, two modifications were made to the model construc- 
tion technique. First, the ice was frozen in a series of 


layers 1-2 inches in depth to minimize the extreme expansion 
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and subsequent fissures. Second, a hollow wooden box 
(Archimedes box) was frozen into the ice immediately behind 
the bollard and at the same nine-inch depth. The Archimedes 
box, which 1S identified here as the core box, reduced the 
weight-to-volume ratio of the model and increased its free- 
board 1.5 inches. The core box and the installation are 
shown in Figures 15 and 16. The combination of layered- 
freezing and flotation box installation resulted in a 
hydrodynamically stable, non-fissured, model with sufficient 
freeboard for testing as Shown in Figure 17. 

The addition of a small bollard on the rear of the core 
box improved the waterborne handling of the model. Further- 
more the box served as a reference centerline and anchoring 
point for test inStrumentation. The only further modifica- 
tion made to the model throughout the remainder of the test 
program was to include the towing bollard as an integral 
part of the core box as shown in Figure 18. The revised 
method of ore the bollard resulted in a weight savings 
of approximately 25 lbs, while reducing the quantity of 
foreign material frozen into that area of the model. From 
all considerations ranging from transportability to hydro- 
dynamic similitude, the model was satisfactory for the 


test series. 


fee GENERAL TEST PROCEDURE 
The test procedure described below is typical of the 


series of experiments. Specific deviations from this 
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general procedure are described in the summary of each 
experiment later in this chapter. 

Approximately two weeks prior to each experiment, the 
ice container was delivered by the Naval Postgraduate School 
Public Works Department to Growers Ice, Salinas, California. 
The container was placed in a freezing tunnel and filled 
with a 1-2 inch layer of demineralized water and allowed 
to freeze. Once the layer was eereee the tunnel was opened 
and allowed to warm to outside ambient temperature, and an 
additional layer of water was added over the frozen ice 
for freezing. This process was continued to the 9-inch 
level at which time the core box, bollard, and thermocouples 
were installed. The process was then continued to the 
18-inch level with the model components firmly frozen in 
position. The average temperature of the tunnel during 
freezing was O°F, and the total time to complete the 
freezing process was 10-12 days. The model was further 
Meeted tO minus 18 -20 °F twenty-four hours prior to pickup 
Mameorder to inhibit melting during the transit and testing 
preparation period as previously discussed. 

In order to exploit the early morning quiesent wind and 
sea conditions, to avoid marine traffic interference, and to 
Minimize the melting by direct sunlight, the tests were 
conducted between 0700 and 1100. The ice was brought to the 
launching site at the U.S. Coast Guard Pier, and launch 


preparations were made. Launch preparations included 
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inspection of ice for fissures, dimensional measurements, 
hookup of instrumentation, and attachment of towing cables 
ama titting rigging as shown in Figures ll, 12, 16, and 18. 
At launch time the entire assembly was weighed and lowered 
into the water as pictured in Figures 19-22. As seawater 
flowed into the container, the assembly was lowered slowly 
by the crane until the buoyant ice model floated freely as 
fe-mein Figure 23. Connections OS made to the towing boat, 
the model pulled clear of the rigging, and the underway 
portion of the experiment commenced as pictured in Figure 24. 
Primarily for safety purposes, considering that the 
tests were to be conducted in a restricted harobr, the 
towing vessel for the first experiment was a Coast Guard 
patrol vessel as shown in Figure 25. All subsequent tows 
were conducted utilizing a Boston Whaler with a small out- 
board motor pictured in Figure 26. A wire bridle connected 
both arms of the towing bollard to a tow rope, and the 
thermocouple wee bundle was streamed parallel to the tow 
rope as illustrated in Figure 27. The ice was towed at an 
average speed of 1 knot (1.78 ft/sec); the speed was 
determined using a Marine Advisor Ducted Current Meter and 
Readout, shown in Figure 28, which was suspended on a rigid 
arm from the towing boat. The model was towed for approxi- 
mately thirty minutes ina "race track" shaped pattern 


approximately 600 yards (.34 mi) in length on a side. The 
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actual time of each tow varied, with the object being to 
have the model underway for a sufficient period to allow 
the full development of surface ablation characteristics, 
while maintaining adequate structural strength for recovery 
at the pier. During the tow, temperatures of the ice 
interior, ice-seawater boundary layer, and model wake were 
recorded. 

Recovery was accomplished by positioning the ice model 
in a web cargo net with a tubular-steel bottom frame and 
raising the entire assembly onto the pier as sequenced in 
Figures 29-31. The assembly was weighed as before to 
determine the total mass lost during the test. Once on the 
pier, the model was photographed thoroughly, surface ablation 
characteristics noted, and overall dimensions recorded. 
Subsequently the ice was broken for disposal and recovery 


of the core box and thermocouples. 


Pee oUMMARTES OF TESTS IN MONTEREY BAY 

The following sections are brief discussions of each 
test conducted in the Bay. The discussion includes details 
of deviations from the general test procedure and evolution 
of the test procedure as experience was gained on each 
experiment. 

1. Test No. 1, September 23, 1977 

The ice was transported from the vendor to the launch 

Site after a four-day freezing period as a single large 


Block. The characteristic "hump" and fissures previously 
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Meeee bed were present as a result of the freezing technique. 
A large wooden reference grid was attached prior to launch 
as seen in Figure 12; the grid provided a reference for 
measuring regression of the ice surface. The container was 
placed in the water, and the ice floated free and was taken 
under tow as pictured in Figure 32. The large towing vessel 
had difficulty maintaining the desired speed of 1 knot. 
The low freeboard and numerous eeures allowed water to 
flow over the top surface resulting in premature breakup of 
the ice and detachment of the towing bollard while underway. 

The bollard and reference grid were recovered; 
however, the ice was not recovered. Therefore, Significant 
data were not obtained. Though this test was largely incon- 
clusive from a scientific view, the photographs obtained and 
the experience gained in the launching and towing procedure 
were useful for future experiments. 

fee test NO. 2, September 26, 1977 

The ice was removed from the freezing compartment 
after a three-day freezing period as a Single layer as before. 
However, three days were apparently insufficient to solidify 
the model completely. Ice formed on the outside with water 
on the inside. The water was covered by a thin layer of ice. 
The thin layer ruptured prior to arrival at the launch site 
as illustrated in Figure 33. The overriding point which is 
apparent from Figure 33 is that the center of the ice did 
not have sufficient time to solidify. The second experiment 


was scheduled only a few days after the first experiment in 
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order to obtain additional information for presentation at 
B@emerrst International Iceberg Utilization Conference on 
Metooer 3, 1977. 

The container with the partially melted ice was 

returned immediately to the vendor to resume freezing. 
feeeerest No. 3, October 3, 1977 

tise bee model for Test No- 93 was frozen during a 
five-day period as a single layer. Note that the ice for 
Test No. 3 was in part the same as for Test No. 2. The 
"hump" and fissures were again present on the top surface. 
The reference grid was attached, and the model was launched. 
The model freeboard was slightly higher as pictured in 
Figure 34 but was still insufficient to prevent water flow 
Over the top surface. The towing vessel used was a 16-foot 
Boston Whaler with an outboard motor which had little diffi- 
culty maintaining the desired speed. Underway photographs 
were obtained; however, the reference grid could not be 
used to determine regression rates over the length of the 
model because of the relatively small differential in these 
rates and the inability to obtain photographs directly over 
the ice. 

The ice model was recovered successfully for the 
mest time as illustrated in Figures 35 and 36. The recovery 
net shown in Figures 35 and 36 did not have the bottom 
tubular steel frame pictured in Figure 30. The tubular 
frame, which was included for later tests, spread the net 


Over the width of the model. Also, the frame submerged 
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the net in the water to ease the positioning of the model 
maeehe net for recovery. In addition, the frame distributed 
the lifting forces over the bottom of the model, thus 
minimizing loads applied to the delicate sides and corners. 

The model was placed on the pier where the weight, 
dimensions, and surface ablation characteristics were 
recorded. Numerous photographs were taken similar to 
Figure 37 which illustrate the erosion caused by the sea- 
water flowing through the fissures on the model. Two signi- 
ficant surface features were noted: First, a heavily eroded 
cavern~shaped area on the bow as shown in Figure 38 and, 
Second, aS illustrated in Figures 39 and 40, a series of 
well-defined, vertical, periodic ripples extending the length 
of either side of the model. These two surface character- 
istics also were observed in subsequent experiments and are 
discussed with other test results in Chapter VI. 

One surface characteristic that was not present on 
this and subsequent tests was the decisive "necking" observed 
on the laboratory models. However, if this process did 
Occur, it was only present in the first few inches and then 
waS SO insSigifniciant that its effect was masked by the 
ripple system on the model. An amplifying discussion of 
this surface characteristic is presented in Chapter VI. 

This test was judged successful in that a complete cycle 
from launch to recovery was accomplished along with obser- 


vations of the ice model dimensions and weight which were 





Measured aS a result of the cycle. However, modifications 
to the freezing technique, increased freeboard, and improve- 
ment of the useless reference grid were necessary for the 
next model test. In addition, temperature sensing devices 
and ripple measuring techniques were required to increase 
the amount of experimental information obtained. 
4. Test No. 4, November 6, 1977 

Test No. 4 marked the atade AO IO) Wjersslilie, Geir 
prehensive data on the ice interior temperatures, ice- 
seawater boundary layer temperatures, and ripple wavelength. 
The model was frozen using the layered technique previously 
described, with the bollard, core box, and thermocouples 
frozen at the nine-inch level. Figure 4lais a top view of 
the model illustrating the placement of these devices with 
Sameamplifying pictorial view in Figure 42. The numerals in 
Figure 4la denote the amplifier channel for that particular 
thermocouple. 

The eee were 1/16 inch diameter Omega 
fast responses Copper-Constantan Type T with beaded sensory 
tips. Five of these thermocouples were bonded together near 
the tips to form an array of five thermocouples spaced at 
1/16 inch intervals (channels 11-15). The array was mounted 
On the bollard (Figure 41b) and placed in the ice with the 
Outermost thermocouple 1.5 inches from the side surface. 
The purpose of the array was to indicate both the interior 
ice temperature profile as the side receded towards the 
array and the ice-seawater boundary layer profile as the 


Meeeoay became exposed to seawater. 
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Interior ice temperatures were recorded on channels 
16-20 (Figure 4la). The output of a particular thermcouple 
was selected by a channel selection and amplified with 
accompanying electronic reference junction (zero degrees 
Centigrade) as shown in Figure 4l1b. The output of the digital 
voltmeter was recorded by voice uSing a tape recorder. The 
Gata were later transcribed and converted to °C using the 
thermocouple conversion tables. . 

The model was launched without incident and, as 
pictured in Figures 43 and 44, had a completely flat top 
Surface and sufficient freeboard. At regular intervals all 
thermocouples were read and recorded. The model was towed 
until all five thermocouples in the array indicated ambient 
Seawater temperature and significant undercutting was observed 
aS in Figure 45. At that point the model was recovered 
without difficulty. Once clear of the water, the extent of 
the differential in melt rates between submerged and above- 
water Bee races readily became apparent as shown in Figures 46 
and 47. Close examination of Figure 47 shows the exposed 
Smeaye pLOtruding down from the top surface. 

Once on the pier, the various layers comprising the 
model were very well defined as shown in Figures 48 and 49, 
as were the ripples and heavily eroded bow section noted in 
the previous experiment. The ripples were measured by 
pushing a multi-fingered contour gage against the surface 
until the fingers conformed to the surface shape. The 


contour was traced onto a strip-chart paper in succession 


47 





with the end result being a profile of pie ice surface 
drawn on the paper. Relative depths were measured using a 
metal reference frame prepositioned on the centerline of 
the core box as illustrated in Figure 50. 

Test No. 4 was considered successful, although the 
ripple measurement method was cumbersome and needed improve- 
ment. The effect of the ice layering on the regression rates 
and ripple formation was a question. In the next test 
the freezing layer thickness was increased to three inches 
to evaluate any dependence on ice layer thickness. [In 
addition, the rapid regression at the ice-seawater boundary 
layer with the subsequent exposure of the array resulted in 
difficulty of monitoring these channels with the current 
Mm=-@ording equipment. For future tests a single thermocouple 
mememeairect output to a portable strip-chart recorder was 
utilized to collect boundary layer temperature profile 
Gata . 

ee Test oe pre December 105-1977 

The alternative 3-inch layer freezing technique was 
applied to this model with adverse results. The increased 
net expansion of the ice buckled the bottom and the sides 
of the core box breaking water tight seals. Broken seals 
allowed water to flow into the box with subsequent freezing. 
As a result the buoyancy of the box was reduced. The top 
of the core box was opened at the launch site, and the ice 
Within the box was chipped out as shown in Figure 5l. The 


box was resealed, and the model was launched. Although the 
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bottom surface of the box was buckled upward reducing volume, 
sufficient buoyancy remained to provide an adequate free- 
board during the test. 

The thermocouple arrangement for the ice model is 
shown in Figure 52. The arrangement was designed to measure 
boundary layer profiles near the bow stagnation point (Channel 
13) and at the indicated side points (Channels 11 and 12). 
Since only one portable strip-chart recorder was available, 
the thermcouples were placed accurately at carefully selected 
distances from their respective receding edges so that they 
would be sequentially exposed and recorded. The sequence 
of channels was 13, 11 and 12. The sequence was determined 
uSing previous regression rate data. The thermocouples were 
monitored initially on a digital readout and then trans- 
ferred to the strip chart recorder when the output approached 
0°C. Each channel was switched back to a digital readout 
when the strip chart indicated ambient seawater temperature. 

immeacdatt1on to the boundary layer, an effort was 
made to measure the temperature profile across the wake of 
the model. The purpose of this effort was to determine the 
extent of cooling experienced by the water flowing past the 
Meee ff cooling proved to be significant, a second ice 
body towed in the wake of a first would be "insulated" 
effectively by cooler water; and the melting rate of the 
second body could be reduced considerably. The temperature 
measurement was accomplished with an array of seven thermo- 


couples positioned across a bridge between two pontoons as 
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shown in Figures 53 and 54. The pontoon bridge was pulled 
behind the ice by a rope attached to the rear bollard on 
the model. Figures 55 and 56 illustrate the towing arrange- 
ment which had the thermocouples positioned one model width 
(4 ft) behind the ice. A distance of 4 feet was a compromise 
based on being sufficiently close to the model to detect 
cooling while remaining clear of the immediate eddies behind 
the model which could distort overall temperature gradients 
in the wake. The array measured the wake temperature at 
a depth of 6 inches, and over a lateral distance 30 inches 
On either side of the centerline. Thermocouple output was 
voice recorded from the digital indicator at periodic 
intervals during the test. 

During the underway portion of the test, the added 
drag of both the pontoon bridge and fouling by floating 
kelp proved beyond the capacity of the electric outboard 
mecor to maintain desired speed. Thus only 80 percent of 
desired speed was attained during the test. The lower speed 
had the effect of decreasing the ablation rate and increasing 
the time required for the boundary layer thermocouples to 
become exposed. Although the boundary layer thermocouples 
eventually became exposed and their outputs recorded, 
recovery of the model was necessary due to a large crack 
forward of the core box precluding exposure of the stagnation 
point thermocouple. 

The recovered ice is pictured in Figure 57 where 


close examination shows the stagnation point thermcouple 
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Gueevtotble through the ice surface, while the forward side 
thermocouple tip clearly protrudes thorugh the side of the 
ice. The crack forward of the core box also can be seen 
clearly in Figure 57. The bow erosion, ripples, and under- 
cutting were again observed in this test and are pictured 
in Figures 57-59. These features were less pronounced, 
however, due to the reduced towing speed. Ripple data were 
obtained using the same contour gage method. An improved 
method was under development but not available for Test No. 
5. The thicker freezing layers also were well defined, but 
their effect on the ripple system and undercutting was 
inconclusive due to the lower speed of Test No. 5. 

A matter of practical importance became apparent 
during this test in that the complexity limitations of the 
ice model testing were being rapidly reached. The towing 
vessel which required an operator, line handler, and data 
recorder/instrumentation operator, plus all the equipment, 
Memermaning Out of room! Also, trying to monitor and record 
Simultaneously, with minimal acceptable accuracy, multiple- 
point thermcouple outputs on two different instruments was 
Serer icult. 

The test waS viewed as a success despite the problems 
with the core box, towing speed, and crack development. A 
larger electric motor was used in future experiments to 
increase the power available to overcome excess drag. The 
freezing technique reverted to the original 1-2 inch layer 


thicknesses. If additional wake temperature data were to 
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be obtained, the pontoon bridge would require modification 
so that it would remain exactly behind the ice model and 
Meemaritt to either side. The pontoon vehicle did not 
follow in the wake of the ice model. In addition, a second 
portable strip chart recorder was obtained to reduce the 
Switching complexity when recording multiple-point boundary 
layer temperature profiles. 
feerest No. 6, April 27, 19738 

In addition to the overall test objectives of previous 
experiments, two additional areas of inveStigation were 
considered during this test. The first area was the 
determination of the effect, if any, of the turbulence 
created by the towing cable and thermcouple wire bundle on 
the regression rate of the bow of the model. The second 
area was the examination of the difference in melt rates at 
two points on a Side ripple. 

The ice was frozen using the 1-2 inch layer thickness 
technique eatin the core box and thermcouples placed as 
Shown in Figure 60. The box had been repaired from the 
damage incurred in Test No. 5, and the towing bollard was 
added aS an integral part of the structure of the box. The 
model appeared in excellent condition at the launch site as 
illustrated in Figure 61, which also shows the thermcouple 
poSitions for this test. The channels were assigned as 
Mmeebows: Stagnation point boundary layer (ll), bottom 
Surface boundry layer (12), ice center forward (13), forward 
ripple boundary layer (14 and 15), rear ripple boundary 


Mayer (16 and 17), and rear ice center (18). 
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Previous ripple data measurements. indicated the 
average peak-to-peak wavelength of the ripples was approxi- 
mately 6 inches. The melting rate varies at different 
points between the crest and trough of a ripple. The arrange- 
ment of thermcouples was an attempt to measure the different 
melt rates due to a ripple. Two stations were selected 
where boundary layer profiles previously had been recorded. 
A pair of thermcouples was positioned at each station with 
a 4-inch separation. Initially the sensory tips were at 
identical distances from the receding edge and at equal 
depths in the ice. Both thermcouples were monitored simul- 
taneously on a real-time basis by a two-channel strip chart 
recorder. Since the separation was less than one ripple 
wave length, the position of each tip relative to the ripple 
was indicated by the sequence of exposure and shape of the 
temperature profile. 

The thermocouples were placed in the ice to achieve 
the following exposure sequence: stagnation point, forward 
Side station, rear side station, and bottom station. Test 
No. 6 was the moSt ambitious attempt at multipoint boundary 
layer recording throughout the series of tests. However, 
Taunching difficulties during the test diminished the thermo- 
couple data. The difficulties began when the ice container 
was placed in the water and prematurely lowered before the 
ice could free itself. This resulted in two adverse conse- 


Guences; first, the thermocouple electrical connections were 


5S 





meeeeea effectively nullifying their output, and, second, 
the ice separated from the container unevenly stern-first 
putting excessive stress on the forward section; as a 
result, the forward section of ice broke at a location 
behind the forward ice center thermocouple as pictured in 
Figures 62 and 63. The decision was made immediately to 
disconnect the thermcouple wires from the forward section 
and continue the test as planned ean the remaining section 
of the model. 

The ice model was pushed from the rear rather than 
towed during this test, as pictured in Figures 64-66, to 
determine the effect of the towing cable and thermcouple 
wire bundle on the development of the bow erosion pattern. 
The wood pushing bar was hinged on either end to allow 
vertical movement of the model and boat. The bar was bolted 
With a facing plate to the rear bollard while clamped to 
the bow of the boat. Two torque lines, running from both 
arms of the main bollard back to the line handler on the 
eae, aided in directional control of the model and mini- 
mized the torque applied on the rear bollard and bow of the 
boat. The underway portion of the test was continued until 
Significant undercutting was observed on the model, at which 
time the ice was recovered as shown in Figures 67 and 68. 

Observations of the model on the pier showed that 
the same bow surface erosion characteristics developed to 
the same extent on this test as when the model was towed. 


Although the bow erosion did not have exactly the identical 
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Shape because of its slanted front as a result of the break, 
close examination showed clearly that the same process of 
erosion had taken place. The well-defined periodic ripples 
were present as in previous tests. An improved method of 
measuring these ripples was utilized as shown in Figures 69- 
71. A wooden reference frame was prefabricated and posi- 
tioned in slots on the core box. A Longfellow Linear Motion 
Position Transducer Model Mecano sicturcd ibigh Ih aegbliex=y M7) (0) 
was mounted on either rail of the frame and moved the entire 
length of the frame is illustrated in Figure 7l. As the 
transducer moved along the rail, the position arm moved 
over the ice surface, and the resulting changes in trans- 
ducer output voltage were recorded by a strip chart recorder. 
The end of the chart paper was attached to the transducer 
chassis so that the transducer moved along the rail at the 
Same speed as the chart paper moved through the recorder. 
Thus the recorder printed an exact geometrical profile of 
the ripples aoe the length of the model. 

This test was the last in the series to study a 
Single ice model. The knowledge and experience gained from 
the series were applied to another test in which two similar 
models were towed in tandem to determine the effect of the 
first model on the melting rate and surface regression 
characteristics of the second. This test and the work 


leading to it are discussed in Chapter V. 


55 





FF. COMMENTS 

An invaluable store of knowledge and experience on 
testing large scale ice models underway in seawater was 
gained as a result of the series of experiments in Monterey 
Bay. While not every test objective was met on every test, 
those that were accomplished added to the experience that 
was nonexistent prior to this series. 

One important point was highlighted by the test series. 
That was the necessity of controlled laboratory conditions 
for experiments with the wide scope of the tests. The 
arrangements alone for the ice model, transportation, launch 
site, boats, equipment, and assistance required the coordina- 
Memenmeor Six different agencies. Add to this the unpredic- 
tability of weather, and one soon grasps the extensiveness 
in preparing for a single test. The complexity involved 
during the tests in monitoring instrumentation while con- 


trolling the boat and model has been previously discussed. 


- 
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V. TANDEM ICE MODEL EXPERIMENTS 


A. BACKGROUND 

The ice model test series described in Chapter IV showed 
that the most intensive melting and subsequent regression 
occurs at the bow of the model. This rapid melting is a 
mgum@etron Of pressure gradients due to variation in the local 
velocity and the higher local temperature gradient profile 
on the bow of the model. The "cavern"-shaped eroded area 
on the lower frontal surface of the model is that area at 
which the higher local velocity and acute pressure gradients 
are most pronounced, which results in the high melting and 
regression rate in this portion of the bow. Griffin's 
[1,2] model showed the strong dependence of melting rate 
on the differential between the ice surface temperature and 
the local ambient seawater temperature. 

One possible solition to minimizing the pressure and 
velocity gradients and the ice surface-seawater temperature 
differential is to tow the iceberg in the wake of a 
"Sacrificial" iceberg. The sacrificial iceberg could shield 
the rear iceberg by lowering the velocity of the seawater 
impinging on the bow and also provide a plume of seawater 
cooled to a lower temperature than the ambient seawater 
temperature. The former would reduce the pressure gradi- 
ents on the bow while the latter would reduce the ice surface- 


Seawater temperature differential, possibly resulting in a 
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Substantial reduction of the high regression rate of the 
rear iceberg. | 

While the general theoretical heat transfer principles 
were well-known and understood, previous research in this 
area with ice bodies underway in seawater was nonexistent. 
As a concluding experiment for this thesis and an introduc- 
tion into possible future work, a test was conducted in 
Monterey Bay consisting of towing two large-scale ice 
models in tandem on December 10, 1978. The models were an 
identical pair and similar in size and construction to those 


in the previous Ssingle-model series. 


eee SORATORY TESTS 

Due to limitations on time and resources, only one large- 
Scale experiment on tandem-towed models could be conducted 
Meretnis thesis. In order to obtain decisive results from 
the large-scale test, a Series of laboratory experiments 
were conducted first on small-scale models utilizing the 
water-channel pictured in Figure 5. The characteristics 
of this water channel are listed in Table I. The primary 
issue to be addressed in these experiments was the distance 
the second model should be towed behind the first to illus- 
trate most drmatically the "protection" the bow of the 
second would derive from the presence of the first. Of 
practical concern was the realization that full-size ice 
bodies could not be towed so close in tandem that longi- 


Pucinal variations in their individual velocities could cause 
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an “accordian" effect which would result in catastrophic 
collisions between the lcebergs. 

Froude number scaling was utilized to model the Monterey 
Bay 1ce model for the laboratory experiments. In Chapter III 


it was stated that for Froudenumber scaling: 
FD nodel D prototype 2) 


while eq. (2) showed that 


with velocity (V) and model length (L) as the controlling 


parameters. Substituting: 
as ey ( =) (10) 
YgL model /GL prototype 
rearranging; 
- Liodel 1/2 
Vmodel ~ Vorototype (F, a) 
: prototype 


igemlength and velocity of the prototype (16 ft., 1.78 ft/sec) 
are known. A freezing container was available with dimen- 
Smons of L = 13.5 inches, W = 4.5 inches, and D = 2.5 inches 
For forming the laboratory model. From eq. (11) the model 


velocity was 0.472 ft/sec. The model Froude number from 
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eq. (2) was 0.023, thus precluding attached bow waves. The 
Bewecnannel velocity eliminated any possibility of signi- 
ficant blockage or wall interference. Although the width 
and depth of the model was not scaled exactly relative to 
the length, their sizes were satisfactory for the laboratory 
tests since drag of the model was not at issue. 

The laboratory models were frozen in a standard freezing 
compartment for a two day period ae a temperature of O°F. 
The model was frozen as a Single 2.5-inch layer with two 
wooden handling devices. Two devices were frozen in each 
Medel tO assure lateral directional stability during the 
tests. 

The laboratory tests were conducted by placing the ice 
models in the flow channel and allowing them to float using 
Pere Own natural buoyancy as pictured in Figure 72. Direc- 
tional control was attained by forcing the wooden handling 
devices to ride ina slot formed by two guide rails that ran 
POngitudinally the length of the channel centerline. Each 
test was conducted for a continuous period of two minutes. 
This two-minute period was determined by placing a single 
model in the channel for a sufficient period to allow 
decisive erosion of the bow as pictured in Figure 73. The 
rear model was placed initially at a separation distance 
corresponding to one model width (W = 1.5 inches), the 
characteristic model dimension chosen for the tests. On 


each successive test the rear model was moved incrementally 
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Buses reducing the separation distance. During each test 
the flow immediately adjacent to the front and rear models 
was Studied by injecting a colored dye in the water in close 


proximity to the area of interest as shown in Figure 73. 


fen eoULTS OF LABORATORY TESTS 

The decisive "necking" observed in the initial laboratory 
experiments was also noted during this series of experiments. 
This "necking" effect served a very useful purpose during 
the tandem laboratory experiments by providing a well-defined 
bow erosion Signature for comparing erosion on the two models. 
Ghapter I and Chapter VI describe the "necking" on the 
laboratory models and the flow characteristics which cause 
Smee rapid erosion to occur. The flow recirculation was 
seen easily with the aid of the injected dye as pictured 
mierigure 73. 

When placed at the initial separation distance of one 
model-width (W), the rear model bow eroded to exactly the 
Same shape as the forward model bow. The rear model appeared 
to have derived no substantial benefit from the presence of 
the forward model. At separation distances of 2W/3 and W/2 
the rear model bow eroded in the identical shape of the 
forward model, but at a slower rate. At a separation dis- 
tance of W/3, the rear model bow eroded in a decisively 
different manner from that of the forward model as shown in 
the sequence of Figures 74-76. The decisiveness of this 


change in erosion pattern can be observed further by examining 
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Figures 77 and 78, which show the bottom surfaces of the 
forward and rear model respectively. The complex flow 
pattern between the two models at the W/3 separation distance 
is shown in Figure 79. 

As a result of this experiment in the laboratory, the 
large-scale tandem-towed model test in Monterey Bay was 
conducted utilizing the W/3 separation distance between 
models. Shorter separation sa SE ARERE were not examined due 
to the practical uncertainties of longitudinal oscillations 


with resulting collisions mentioned earlier in this chapter. 


SeeerMONrTEREY BAY TANDEM TEST, DECEMBER 10, 1978 

The models used in the tandem experiment were frozen 
in identical fashion, utilizing the 1-2 inch layer freezing 
technique with similar core boxes installed. Thermocouples 
were not installed in the ice due to the complexitv of the 
experiments in handling two large ice models simultaneously. 
Since only one crane was available, the models were launched 
separately (forward model first) and then linked together 
in the Bay as pictured in Figure 80. The linking was accom- 
plished by two rigid metal rods connected to either arm 
of the towing bollards that were sized to provide a constant 
W/3 (16 inch) separation distance between the models. 

The models were towed over the same course aS in previous 
experiments utilizing the same Boston Whaler towing vessel 
as pictured in Figure 8l. While the rear model remained 
behind the lead model for some brief periods (Figure 82), 


lateral instabilities in the models resulted in lateral 
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motions and misalignment of the models during the major 
Merton of the experiment as shown in Figures 83 and 84. 
Attempts at restoring the misalignment by providing a drag 
force from a following Boston Whaler on the model train 
(Figure 84) were unsuccessful. 

Recovery was initiated after 21 minutes underway when 
Significant undercutting was observed on the lead ice model. 
The rear model was recovered first, followed by the lead 
model. Both models were again weighed to determine total 


net weight lost during the test. 


feeeeeeoULTS OF TANDEM BAY TEST 

While decisive indications were obtained from the 
experiment, the results were ina large measure inconclu- 
Sive due to the nature of the experiment. Weight indications 
Showed that the lead model lost 2750 lbs. and the rear model 
2000 lbs. While this large weight loss differential would 
appear to prove decisively the "sacrificial" iceberg theory, 
One must remember that the lead model was in the water an 
additional 36 minutes due to the limited single launch and 
recovery capability available. 

The recovered lead and rear models are shown in their 
bow view in Figures 85 and 86 respectively. The bow erosion 
of the rear model is approximately the same as the lead 
model which would negate the shielding theory provided by 
the lead; however, the rear model was underway eight addi- 


tional minutes over the lead in order to link up with the 
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Meeutally drifting lead model in the Bay. Additionally, 
the misalignment resulted in the rear model's bow being 
exposed to the full flow of the oncoming seawater during the 
major portion of the test. 

The Tandem Bay test procedure can be improved by 
utilizing two lifting cranes to effect a simultaneous 
launch of both models and by connecting the tandem models 
through a lightweight rigid frame 55 prevent lateral motion 
between the models. The former would be more expensive and 
require additional pier space, but the result would be a 
more accurate reflection of the actual underway mass loss 
differential between the models. The latter would assist 
in better defining the bow erosion differential between the 


models. 
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Vine EXPERIMENTAL RESULTS WITH THEORETICAL CORRELATIONS 


A. GENERAL 

The following are discussions of the experimental 
results determined from the ice model test series described 
in Chapters III-VI. Comparisons with appropriate theoretical 
predictions or other experimental results are made including 
descriptions of the theoretical models and experimental 
techniques. 
SeeeecOMPARISON OF EXPERIMENTALLY MEASURED ABLATION RATES 

heed THEORETICAL PREDICTIONS 

Griffin developed a model for predicting the ablation 
rates of icebergs moving through seawater in laminar flow 
eeeeeHe modified his model [2] to account for turbulent 
flow in conjunction with the large-scale ice model experiments 
Bemauected at NPS. 

The analysis model considers turbulent flow over a 
flat plate with zero pressure gradient and solves the boun- 
dary layer equations for simultaneous heat, mass, and 
momentum transfer. The two-dimensional problem is reduced 
to one dimension by integrating across the boundary layer. 
Mies phase transformation from solid to liquid occurs under 
Steady-state conditions enabling a coordinate system to be 
affixed to the melt interface. A constant density is 
assumed, and the ice is a pure solid phase. Gravity is 


neglected, thus neglecting the buoyancy effects on the heat 
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Meeeretr. The model requires inputs for the physical 
Memerties Of Seawater and ice, including the applicable 
Semensions and free-stream velocity. The properties of 
mee and seawater are assumed to be constant for the solution. 

The graphical results for overall surface regression 
for three of the large-scale experiments in Monterey Bay 
are shown in Figures 87-89 and summarized in Table V. 
A comparison of the actual regression rates at three points 
Memeene model tested on October 2, 1977, with the prediction 
from Griffin's Computer Code is shown in Figure 90. The 
agreement between experimental values and theoretical pre- 
dictions is good. A plot of turbulent flow Stanton numbers 
versus Reynolds numbers is shown in Figure 91 for the 
range of Prandtl numbers typical for seawater. Data points 
from Table V of NPS experiments are shown again to be in 
good agreement with the theoretical curves. A similar 
plot of Stanton numbers versus Prandtl numbers is shown in 
Figure 92 oso range of Reynolds numbers from laminar to 
meeoulent flow conditions. 

The conclusion drawn from the above comparisons was 
that the Griffin turbulent-flow iceberg ablation model was 
Satisfactory for predicting the ablation rate at various 
points on the sides an iceberg given the appropriate physi- 
cal conditions such as ice and water properties, velocity, 
and water temperature. However, the bow area of the iceberg 


requires further treatment. 
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C. MEASUREMENT AND VERIFICATION OF THERMAL BOUNDARY 
LAYER PROFILE | 


The profiles of the water temperature versus distance 
from the ablating surface for experiment Nos. 4 and 5 in 
Monterey Bay are shown in Figures 93-95. The data on 
Experiment No. 4 were obtained by recording a digital readout 
from the thermocouple array described in Chapter IV while 
the data for Experiment No. 5 were recorded by a single 
thermocouple output to a continuous strip chart recorder. 

Due to the integral solution technique applied across 
the mass diffusion, thermal, and momentum boundary layers, 
Griffin's model does not provide predictions for the 
thermal boundary layer profile. However, an approximate 
determination of the accuracy of the initial profile slope 
could be calculated by applying a one-dimensional energy 
balance across the ice-water interface as shown in Figure 
O07 . 


From Figure 107: 


P. ean : 
qd, = heat transfer from water (opens 
= heat transfer into solid ice (he 
ea Seccm 
M =~eelatent neat Of fusion (eae) 
The energy balance is: 
or eer ae Wi ae qs (BE), 
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where 


oa ice density Con? 
Cm 
Ve = regression rate (=a 
Also: 
Shy Se ee a) (13) 
and 
ae Sy Ks, = ; (14) 
where K is the respective thermal conductivity and = is the 


respective initial slope of the temperature profile. 


fMeem Ref. [7] at 00°C: 


R = 6.35 8 on acca 
at Crem sec 
= ame cal 
Sy aS) See ae ae SGnGnacaa 
Mie =e 91977 cee 
g 
2 
9 = .9168 g/cem 
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The ice interior temperature Profile - Test No. 4: 


From Figure 96; 


ar 
oy 


Wl Oe. 
. Mees an) (2.54 2) 
ale @ 


Cc) 


From Figure 94; 


0, en 
a BOL tn) @o5n Ss 
in 
From Table V; 
V —— ye omlenivadte (at 15 = 2.4 Ia) 


Berving eq. (14) and substituting into eq. (12) gives 


a 
sec-cm 
Solving eq. (13) gives a value of 
— 0.3 ean 
sec-cm 


Reasonable agreement occurs between theoretical and observed 
results indicating that the measured profile 1s a reasonably 
accurate representation in at least the critical initial 


Slope. 
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See RoPPLE EFFECTS AND MEASUREMENT 
1. General 

Throughout the test series, ripples were observed 
to have formed on the sides of the ice models tested in 
Monterey Bay. The ripple formation 1S not unusual as turbu- 
lent flows in nature are observed to produce several differ- 
ent types of regular waves, including wind-generated water 
waves. Carey [8,9], Larsen [10], Ashton and Kennedy [11,12], 
and Hsu [13] have described and analyzed the occurrence and 
behavior of regular waves on the underside of river ice 
covers. These waves are known as ice ripples and can be 
expected to form at any ice surface exposed to a turbulent 
water flow, at least when the boundary is receding due to 
melting. Tatinclaux and Kennedy [14] developed an analyti- 
cal model for formation of ripples on icebergs. 

meee Ripple Effects on Drag 

Both heat and momentum are transferred between the 
ice-water boundary layer and the surrounding flow by means 
of turbulent exchange processes. In addition, the ripples 
eventually will become sufficiently high so that the flow 
Will separate near the crest of each ripple. The normal 
relatively low pressure at the crest will then continue 
across the rear face of the ripple until the flow reattaches 
at the frontal face of the next ripple. The large pressure 
differential across the ripple results in a significant 


drag force exerted on each ripple. Also there will be 
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gravity waves generated in the water by the ripples, thus 
adding to the iceberg drag. 

several methods have been proposed for calculating 
the composite friction coefficient by Carey [8], Larsen 
[10], and Uzuner , et. al. [15]. Larsen showed that drag 
increases almost linearly with ripple steepness (peak-to- 
peak ripple height divided by ripple length) up to a value 
@eeabout 0.1 and much more rapidly thereafter. For fo ty 
developed ice ripples, values of steepness in a range of 
0.-0.15 are not uncommon. Further, such ripples could 
increase frictional resistance by 50-100 percent. 

@eestpple Effect on Heat Transfer 

The formation of ripples at the ice boundary has 
the effect of increasing the rate of heat transfer at the 
boundary. The increase in heat transfer is due, first, to 
the increased hydraulic roughness which affords a larger 
net area available for heat transfer and, second, to the 
increased turbulent intensity and flow separation caused by 
the ripples. 

4. Measurement of Ripples on NPS Tests 

Chapter IV described the various tehcniques used to 
measure the ripples on the sides of the ice models tested 
in Monterey Bay. From the data gathered, histograms of 
ice ripple wavelength were constructed and are shown in 
Figures 97-99. From the histograms the average wavelength 


was determined for each experiment. 
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Experimental programs conducted at the Iowa Insti- 
tute of Hydraulic Research on the formation of ice ripples 
have been reported - Ashton and Kennedy [12] and Hsu [13]. 
The experiments were conducted on an ice slab of dimensions 
me 40 ft. and W = 2 ft. over a velocity range of 0.30 ft/sec 
to 2.7 ft/sec. Hsu reported that the average wavelength 
(A) 1S inversely proportional to the mean flow velocity 
(U) and that the Reynolds number based on wavelength is 


Senstant and given by: 
R = ne = Be 2 2x 10 (aS) 


The standard deviation of the experimental results from this 
average is 1l.l x LO. 

The average wavelength for each NPS experiment was 
determined from eq. (15) and compared to the corresponding 
average wavelength determined from the histograms. The 
results are listed in Table VI and show good agreement. 

See sohORETICAL PREDICTIONS AND MEASUREMENT OF ICE 

elo lLOR TEMPERATURE PROFILES 

The temperature profile within the ice model was 
measured during the underway portion of Test No. 4 in 
Monterey Bay. The thermocouple array discussed in Chapter IV 
and shown in Figures 41 and 42, was read at specific time 


intervals for a period in excess of 37 minutes; readings 


were volce recorded. 
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TABLE VI 


COMPARISON OF RIPPLE DATA MEASUREMENT 


DATE: unOV 7 77 27 EO 97, AG OMIT Fit 
Melocity (U) 5 ies 1.9 
(ft/sec) 
Wavelength 5.4 Get 4.4 
(A,) (in) 
Wavelength 5.4 Seis: el 
(A) (in) 
Note: 


r Wavelength from eq. (15) 


ho Wavelength from histograms. 


74 





Mae transient thermal conduction equation for one 
Spatial dimension was solved, and the results were used to 
correlate the experimental meaSurements. The equation and 
development of the solution is given in Appendix B. Included 
in Appendix B is a Hewlett-Packard-9831-A computer program 
to generate the theoretical temperature profiles, thermal 
Metrusivity of ice (a), and ice regression rate (Vi). 

The comparison of the theoretical and experimental 
temperature profiles is shown in Figures 100-105 for each 
respective time interval. The theoretical curves represent 
a best parametric fit with the experimental values by 
independently varying such parameters as thermal diffu- 
Sivity and regression rate as described in Appendix B. The 
resulting theoretical thermal diffusivity and ice regression 
rate from the best fit condition were 3.0 x ome aC SRE and 
m2 xX to - m/sec respectively. The actual thermal diffu- 
Sivity of ice determined from Ref. [7] was 5.3 x Loe m-/sec 
for pure solid ice. The actual regression rate, determined 
by initial and final ice model dimensions for the 37 minutes 


2 m/sec. 


mee underway time, was 3.53 x 10) 
Figures 100-105 show compatability between the theoreti- 
cal and experimental data for temperature profiles within the 
ice. The theoretical curves and experimental values do 
diverge to some extent at time intervals 5 and 6. One 
hypothesis for the divergence is the significant change in 


thermal diffusivity in the relatively steep temperature 


mradient. 


TS 





The variance between theoretical and actual ice thermal 
diffusivity 1s explained by Mellor [16] as the result of 
air bubble entrapped in the ice which lowers the thermal 
diffusivity from that of pure bubble-free ice. There was 
@ Significant bubble population throughout each model tested 
mn the bay. 

Since the time over which the regression rate was calcu- 
lated did not include the 7 aimates ne ice was motionless 
in the water, it was not surprising that the measured 


regression rate was higher than the theoretical prediction. 


FP. MEASUREMENT OF ICE MODEL WAKE TEMPERATURE 

The method used to measure the wake plume on Test No. 5 
in Monterey Bay is described in Chapter IV. The method 
proved largely unsuccessful for the following reasons. 
Mest, the pontoon bridge did not remain aligned continu- 
ously with the centerline of the model throughout the test 
which limited the opportunities to obtain data. Second, the 
thermocouples had to be read individually in stepwise fashion 
due to recording equipment limitations thus adding variance 
to the recorded temperature profile at any given time. 
Third, variations in surface water temperatures, large 
scale eddies in the wake, wave motion and similar effects 
added to scatter in the measured values. The recorded 
temperature profiles are shown in Figure 106 for five 
different time intervals. A search of the literature showed 


no data for similar measurements for comparison with the 
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NPS experimental results. Further, efforts at measuring 
Mane wake plume were not attempted due to Tie time and 
resources. 

The conclusion is that measurement of wake temperatures 
is not as straightforward as one might expect. While the 
equipment recording problem could be solved by continuous 
monitoring and recording techniques, the large scale eddies, 
which were observed in the wake with dye injection, signi- 
ficantly complicate the interpretation of temperature 
profiles obtained. 

One obvious reason for measuring the ice model wake 
temperatures is to aid inthe determination of the effective 
cold-water "insulation" provided by the model for a similar 
model being towed in the wake. 

One alternative method for measuring the wake tempera- 
ture would be to measure it with the second model in tow. 
The close proximity of the bow of the trailing model to 
the stern of the lead model would have a stabilizing effect 
on the flow leaving the lead model. An attempt at this 
method was planned for the tandem test in Monterey Bay but 
could not be executed because of the lateral oscillations 


of the models. 


Se FLOW SEPARATION ON THE ICE MODELS 
The decisive "necking" observed on the small-scale 
laboratory ice models was due to the recirculating flow in 


the forward one-third region on the models. The recirculation 
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was a reSult of the initial flow separation around the 

sharp frontal corners and subsequent reattachment as des- 
cribed in Chapter III. The separation region, where the 
necking occurred, was observed visually with use of dyes 

as described in Chapter V. However, the "necking" was never 
observed on the large-scale ice models tested in Monterey 
Bay, despite the fact that flow separation around the initial 
90-degree corners certainly aseureaady: Meauscussilon Of. che 
flow separation conditions for the two modeling situations 
Follows. 

Consider first the frontal corners on the small- and 
large-scale models. The corners are initially areas of 
intense heat transfer and erode rapidly once the model starts 
moving in the fluid medium. As the corners erode, their 
radius of curvature increases thus decreasing the tendency 
for flow separation. The actual rate of decrease in the 
flow separation tendency will vary depending on the laminar 
Or turbulent sondelen Greener leva on the Laboratory Cests, 
the flow was laminar due to the low Reynolds numbers involved. 
The small frontal area of the laboratory model did not offer 
Merrclent opportunity for the flow to reach turbulent con- 
ditions prior to rounding the corners. Laminar flow is 
more sensitive to flow separation even on relatively stream- 
lined bodies. Therefore, although the corners became 
rounded on the small-scale model, the flow, which remained 


laminar, continued to separate with the subsequent "necking" 
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as previously seen in Figures 6 and 7 of Chapter III and 
Pagure 73 of Chapter V. 

On the large-scale model, the flow was turbulent as 
evidenced by the full-length ripple system on the sides and 
"Ccavern'-Shaped area on the bow of the model. Due to the 
Significantly larger frontal area on the large-scale model 
as compared to the small-scale model, the flow has an 
enhanced opportunity to reach a turbulent condition prior 
'to rounding the corners. A well-known fact is that turbu- 
lent flow has a lower tendency than laminar flow for separa- 
tion. Therefore, one hypothesis is that once the corners 
on the large-scale model become rounded, the flow separation 
ceased and no observable "necking" occurred. Such a hypothesis 
could be tested using both potential flow theory and actual 
continuous observation of the large-scale model as was done 
in the laboratory experiments. 

The significant difference in the melting shapes of the 
small- and large-scale models emphasizes the importance of 
careful treatment of the scaling factors such as Reynolds 
number when modeling a large iceberg. The melting shapes 
were obviously strongly dependent on Reynolds number, which 
Was on the order of 10° in the laboratory tests and Oe” 


in the Bay tests. 


79 





) tie RE CCUIMENDALEONS FOR FUBURE WORK 


A. GENERAL 

There are many areas of research concerning glacial 
ice moving through seawater which strongly merit future 
examination. Included are improvements in testing techniques 
and instrumentation used to obtain such basic data as drag 
force, velocity, temperatures, and Salinity. The following 


are brief discussions concerning these research areas. 


eeroOriCS FOR FUTURE RESEARCH 
meee asurement of Drag 

BOr a Square bow, such as shown in Figures 27, 29, 
and 31, the drag is dominated by form drag. At the low 
Froude numbers, the wave drag is negligible. Skin friction 
is relatively insignificant. 

As the iceberg melts, the shape tends to a geometry 
yielding lower form drag. Also the bow may be shaped to 
give lower form drag. In that case, skin friction becomes 
More Significant. There iS an obvious interaction between 
Convective heat transfer and skin friction. One expression 
Se the relation is Reynold's analogy. 

As a result of the past towing experience, the NPS 
team should develop techniques to measure drag accurately. 
Drag information is useful in the interpretaion of shape 


Changes and for correlation with melting rates. 
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2. Tandem Towing 

Only one test was conducted at NPS in this inter- 
esting and complex method of moving icebergs. One potential 
advantage of tandem towing is the fluid mechanical insulation 
(cold water) provided by the lead iceberg. The lead iceberg 
can be considered to be a "Sacrificial" iceberg providing 
the cold water for insulation. Evaluation should include 
boundary layer temperature profile measurement and flow 
field surveys at the rear of the lead iceberg. 

fume uence Of Buoyancy on Heat Transfer 

The paper by Josperger [17] clearly demonstrates 
the importance of buoyancy forces on natural convective 
heat transfer. The critical Grashof number for transition to 
turbulent flow is a factor of 10 lower for fresh water ice 
melting in seawater than for a vertical hot plate in air. 

An analytical model should be developed to determine 
the influence of buoyancy. A bucyancy term should be 
introduced in one equations of motion developed by Griffin. 
The impact of this term on heat transfer should be examined. 

The Naval Postgraduate School is about 60 miles from 
San Antonio Reservoir which is near King City, California. 
In order to test the theory concerning buoyancy and heat 
transfer under forced convection (1.e., moving iceberg), 
several icebergs could be towed in San Antonio Reservoir. 

A comparison could be made between melting rates in fresh 
Water and Seawater. Furthermore, the inverse melting 


problem, i.e., an iceberg of frozen seawater melting in 
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fresh water, can be studied to gain insight on buoyancy 
effects in forced convection. 


mesnstabllities of Insulating Material, e.g., 


Flagwaving Instability 


Consider an insulation scheme using thin plastic 
sheets. A protective layer of cold water develops between 
the iceberg and the plastic sheet. 

There are several Are SGN EM nooks, eS Oeises iat [1 Os, 
Chandrasekhar [19], and Shen [20], which discuss stability 
of flows. In rotating coordinates, the book by Greenspan 
[21] is especially relevant. One well-known instability 
Mmenat which occurs in parallel flows of fluids with 
unequal densities. 

A survey is needed of the instabilities that may 
arise in the context of icebergs with insulation. Such 
a survey 1S proposed here. 

5. Solid Mechanics of Iceberds 

The Naval Postgraduate School has a Chair in Arctic 
Meee. The incumbent starting in July, 1978, is Dr. 
Wilford F. Weeks, who is a Research Glaciologist with the 
U.S. Army Cold Region Research and Engineering Laboratory. 
Dr. Weeks has agreed to devote some of his time to the 
question of solid mechanics of icebergs. 

He lS apparent from the technical papers at the 
Iceberg Utilization Conference that many of the existing 
engineering techniques relative to cracks, fractures, stress 
distribution, multiaxial stresses, etc., are applicable to 


iceberg structures. 





6. Development of Instrumentation 

To verify the computer code of Griffin, measurement 
of temperature, salinity, and velocity profiles would be 
especially helpful. 

Velocity of the iceberg is measured uSing standard 
oceanographic equipment. Drag is measured also using 
standard equipment. | 

Measurement of salinity and velocity profiles 
requires special consideration. The surface of the ice is 
receding. The iceberg is a remote vehicle. 

As a research task, techniques should be developed, 
or improved, for measuring: 

Mm velocity of iceberg 

me drag 

eee salinity profile 

ee velocity profile 

5. temperature profile 
As a further aid to the experiment, underwater photographs 
and underwater movies of an iceberg in motion should be 
made. 

[ee Gaboratory Tests 

Towing tests in Monterey Bay are time consuming and 
expensive. Hence, each test must be planned carefully and 
Overall test objectives must be defined clearly. To supple- 
ment the towing tests in Monterey Bay, tests can be conducted 
using the water tunnel and water channel located in the 


Department of Mechanical Engineering. 





By uSing a larger test facility, larger test 
Specimens, and higher flow velocity, turbulent flow may 
be achieved in the laboratory. Special emphasis should be 
placed on defining the usefulness of boundary layer trips. 

Due to the relative simplicity of laboratory testing, 
Many tests can be conducted in the laboratory compared to 
larger scale towing tests in the ocean. 

8. Field Tests 

There have been three options for large scale tests 
discussed and evaluated by the Naval Postgraduate School 
research team. Option 1, which is towing in Monterey Bay, 
has been exercised six times. Option 2 involves use of 
tidal currents at Moss Landing. Elkhorn Slough has a 
Narrow Opening; the bridge for State Highway 1 crosses the 
Memon at this point. Due to tides, a velocity of 1.0 knot 
Becurs and persists for about 30 minutes. Using the bridge 
aS a point to tie the two ropes, tests could be conducted 
with the 16-£00t iceberg. The bridge and the nearby 
navigational pilings could be used as observation points. 

Option 3 is to tow the large scale iceberg in San 
Antonio Reservoir. The goal of these tests would be to 
determine influence of buoyancy on heat transfer. Also at 
San Antonio Reservoir water temperature would allow use of 
Swimmers near the iceberg; the swimmers could observe the 
flow at close hand. 

Experience has shown that the field tests would be 


facilitated greatly if an observer could ride the iceberg. 
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Hence a double wall box 1S proposed. For the size of manned 
iceberg, the weight of the ice is 7735 lbs. With a payload 
of 1500 lbs, which includes the weight of the inner box, 
the freeboard is almost 12 inches. The design uses a broad 
beam to avoid stability problems if ice were to break away. 

Advantages of a manned iceberg include the following: 

ime abllity to observe the flow. 

mee ability to conduct eee ence, start recorders, 
and vary profile instruments. 

9. Insulation 

Many different forms of insulation have been proposed 
including plastic sheets, frozen wood chips, plastic foam, 
etc. An alternate scheme might be the use of air bubbles 
at locations with high heat transfer rates. 

To examine the concept, an experimental program 1s 
proposed. First, the bubble insulation would be investigated 
in the water channel. Second, one or more of the large 
scale Meebergs would be dedicated to an evaluation of 
bubbles for insulation. Third, depending on the outcome 
of the experiments, a study would be made of pumping require- 
ments. Also, the study would examine installation of 
bubble source on the iceberg. 

10. Potential Flow Computer Code 
| Pressure gradients have a pronounced effect on boun- 
dary layers. A potential flow computer code would be used 


to calculate pressure coefficients and pressure gradients. 
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ine Shape of the body will be that shape observed during 
experiments conducted in Monterey Bay. 

A computer code has been developed by Professor 

| C.J. Garrison under a NSF grant. Professor Garrison is a 
faculty member in the Mechanical Engineering Department at 


NPS. 
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AEPENDIX A 


AREAS OF CONCERN WITH ICEBERG DETERIORATION 


iaeetollowing 1S a discussion by Stolfi, et al [22] of 
the motivation and scope of interest in iceberg deterioration. 

A large percentage of the land surface of the globe 
consists of arid and semiarid regions which support less 
agriculture and smaller populations than other areas. [In 
certain semiarid regions where agriculture is successful, 
that activity is dependent upon well-developed irrigation 
technology and the light variable rainfall of the region. 

The arid regions, and especially the areas in them lying 

far from the world's oceans, present extraordinarily challenging 
problems for development in some future time frame. Certain 
arid regions and many Semiarid regions which lie close to 

the world's oceans contain Significant populations and 

present special opportunities for expanded development and 
insurance against the periodical, economically damaging 

droughts associated with rainfall scarce areas in which 
Productive agricultural activity is being conducted. 

Several appraoches can be envisaged for the support of 
agriculture and associated populations and industry in arid 
and Semiarid regions. The light, existing rainfall can be 
gathered, stored, managed, and metered out to the surrounding 
countrySide and cities inwater storage and irrigation works 


which conserve and distribute the existing rainfall. Similar 
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water storage and distribution works can transport water 
from heavier rainfall regions located economically close 

of the drier areas, to make the water available for agricul- 
ture and maintain the water tables and other underground 
water resources. In cases of deep and extensive drought, 
however, aS in Western North America, 1975-1977, the lack of 
rainfall will empty reservoirs, lower permanent rivers, and 
drive down water tables over large regions threatening 
immense damage to grazing, agriculture, and ultimately 
industry. In both arid and semiarid regions in which 
productive agriculture has been developed through water 
conservation, storage, and distribution technology, the 
water supplies can be supplemented by (1) inter-region trans- 
fer of water from distanct heavier rainfall areas, (2) 
fossil-fuel, nuclear, or solar desalinization of ocean 
water, and (3) transfer and terminal processing of Antarctic 


ice. The approaches can be summarized as follows: 


od 


Support OF Agriculture and Industry 


in Water Deficit Regions 











Arid Regions Semlarid Regions 
(Sahara, Namib, Atacama, etc.) (California, Argentina, South Africa, etc.) 
Transfer of Water from Closely Conservation, Storage, Distribution 
Adjacent Areas with Existing Rainfall of Light Rainfall Resources 
Inter-Region Transfer Ocean Desalinization Antarctic Ice 
Of Water from Distant (Using Fossil, Nuclear Ocean Transport 
Areas or Solar Energy) and Terminal Processing 
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Mie expansion of economic activity into arid regions and 
the intensification’ and assurance of uninterrupted activity 
in semiarid regions demand advances in technology related 
with Inter-Region Transfer, Ocean DeSalinization, and 
Antarctic Ice Transport and eventual assignment of financial 
and economic resources to such projects. Inter-region 
transfer involving the movement of water from assured, 
excess Supplies, for example, in the Northwest Territories 
in Canada to the California-Nevada-Arizona area in the 
United States, would involve immense cost, formidable 
engineering difficulties and long time periods. Ocean 
desalinization, especially with the use of nuclear or solar 
technology, offers substantial opportunities for support 
of economic activity although at immense expense, with 
potential physical dangers and extended time of development. 
The transport of Antarctic ice and the melting of it in 
mera and semiarid coastal regions, although a relatively 
new idea and largely untried, offers possibilities for the 
Mapport Of economic activity. The transport of Antarctic 
ice presently appears to be comparatively inexpensive, 
involves the use of existing technology, and is capable of 
implementation in the near future for areas south and possibly 
Seven mOrth of the equator. 

In arid regions, e.g., Sahara, Saudi, Atacama, Peruvian, 
Kalahari, and Australian Deserts, most parts of which lie 
far from adequate rainfall areas and desert rivers, the 


expansion of civilization will be dependent ultimately on 
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Ocean deSsalinization and the transfer of Antarctic ice. In 
semiarid regions, e:g., South Africa, angola, Argentina, 
Iran, Pakistan, California, Mexico, Australia, etc., economic 
expansion and the support of larger populations will depend 
On exploitation of the local rainfall in combination with 
ocean desalinization, the transfer of Antarctic ice, and the 
expanded use of adjacent but distant stored rainfall. The 
western half of the United States, except for the relatively 
small, rain-rich Pacific Northwest stands as a monument to 
agricultural opportunity which is dependent upon the availa- 
bility of water. The impoSing total of 15 western States 
is predominately semiarid or arid and dependent upon irri- 
gation works and scanty rainfall for the sometimes intense 
but mostly modest or nonexistent development of agriculture.* 
Additional immense regions of the globe in similar favorable 
latitudes lie agriculturally dormant because of the water 
deficit and high evaporation and transpiration character- 
istics of the areas. 

In arid and semiarid areas lying south of the equator 
Meee within reasonable proximity to the oceans, e.g., parts of 
me Republic of South Africa, scientific theory and experi- 
ment Support a view that the transfer of Antarctic ice is 


limited more by the problems of potential resettlement and 


eizona, California, Colorado, Idaho, Montana, Nevada, 
New Mexico, North Dakota, Oklahoma, Oregon, South Dakota, 
Texas, Utah, Washington, Wyoming. 





Bme cost of irrigation works than by limits on the ice 
transfer technology. Given the relatively short distances 
from the tabular icebergs, and the comparatively cool 
intervening water and air temperatures for movement, the 
NPS Scientists feel that Scientific feasibility is close to 
being established for movement from Antarctica to South 
Africa, Argentina, Atacama, and Australia. Opportunities 
exist for the United States gavemmentt to increase its 
influence peacefully through the use of Superior technology 
in Southern Africa, South America, and Australia. Engineering, 
economic, and political studies will be required to guide 
U.S. government policy makers in the exploitation of the 
Scientific opportunities for the settlement of populations 
and the -expansion of agriculture and industry in previously 
excluded arid regions. 

The stages of transfer and use of Antarctic ice run in 
the following generalized pattern: 


AnNGarcere dice: 
stages of Transfer and Use 


South ———_—_— ee North 


Locate Breakout Transport Transport Terminal (Extract) 
Select Hookup Insulate Terminal (Dock) Distribution 


The stages through which an Antarctic ice transfer 
project must necessarily proceed are common for ice delivered 


mOvany location on the globe. The individual conditions of 


97 





@ecan routing, terminal docking, etc., may vary enormously, 
for example, between ice delivered to Southern California 
and South Africa, but the general process, or stages, of 
transfer will be virtually identical. With two noteworthy 
,exceptions, the technology for each of the stages is largely 
in hand. Space photographic technologies are available to 
locate and select tabular bergs which have the sizes and 
shapes to meet other delivery criteria. The breakout of 
large glacial tabular bergs may not be necessary depending 
upon the time of the year, e.g., Antarctic summer, and the 
location of the glaciers relative to open water. The hookup 
and towing of the tabular bergs can be forecast to be 
engineering challenges of the first magnitude but surmountable 
with present technology. Ocean routing and terminal environ- 
mental impact should be able to be mastered with existing 
technology. 

The most imposing scientific challenge 1s the one of the 
melting of ice moving in seawater, and little theoretical or 
experimental knowledge is available on the subject. If 
rates of melt can be described adequately and predicted, the 
Scientists can present the ultimate reality of ice transfer 
projects in terms of (1) the amount of melting and (2) the 
associated necessity of lack thereof for insulation. The 
MOst imposing potential problem in the transfer of Antarctic 
lce can be forecast to be the one of placing insulation on 
meme ice and maintaining it there during the transportation 


to the terminal area. If the experimentally supported 
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theory shows that insulation is not necessary for the move- 
ment of tabular bergs to areas like Southern Australia, 
Argentina, and Southwestern Africa, the potential costs of 
movement would fall steeply, and the economics of transfer 
would be more attractive. The absence of insulation reduces 
costs and engineering difficulties to a point where the NPS 
Scientists tentatively conclude that the transfer of ice to 
Australia, Argentina, and South Africa is technically and 
economically within the grasp of governments and the 

larger private businesses associated with those areas. The 
perception of the necessity for agricultural expansion on 
the parts of the governments and the economic atractiveness 
on the part of private business remain to complete a picture 
of Antarctic ice capable of being transferred to southern 
regions in the immediate (one to three years) future. 

Mire transfer of Antarctic ice to regions north of the 
ecuator present greater cost and engineering problems largely 
because of the problems of melting and the prima facie 
necessity for insulation. No realistic schemes presently 
Sexist for the insulation of the large tabular bergs in terms 
of physical or chemical barriers between lice and seawater. 
The NPS scientists, based on recognition of the fact that 
increased water temperature exponentially increases the rate 
of melt of ice, are gathering experimental data both in 
Monterey Bay and in the NPS mechanical engineering laboratories 


On the cooling effects of ice moving through water. A 
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"olume" of cool water is assumed to exist behind ice in 

water which could reduce water temperatures to the degree 
that other ice bathed in the plume would experience signi- 
ficantly reduced melting. The potential insulation scheme, 
in effect, would involve either (1) towing several bergs 
maecolumn with a small "Sacrificial" berg in front or 

(2) towing several bergs in column while removing the lightly 
compacted snow and ice.from the top of the leading berg and 
depositing it in front of the leading berg as a coolant for 
the most exposed ice. 

Several challenging but surmountable problems exist at 
the terminal area for the tabular bergs. The ice upon 
arrival will have imposing drafts of approximately 200 - 500 
feet depending upon the distance and velocity of the tow 
and the water temperatures encountered. The immense drafts 
of the bergs will demand carefully selected terminal areas 
in order to place the ice conveniently close to the receiving 
water system and reduce the impact of adverse sea conditions 
on the floating reservoirs of fresh water. Ice, for example, 
delivered to Southern California could be "docked" immediately 
memeneor Santa Cruz Island, 20 miles from the California 
coast protected from the prevailing northwesterly winds and 
representing a reduced impact on the local environment. The 
most difficult engineering problem at the terminal location 
would be the one of establishing a barrier between the tabular 
bergs and the seawater in order to extract the most fresh 


water possible from the ice. Given the static situation of 
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the bergs at the terminal location, one can envision 
engineering solutions Banging from the production and trans- 
fer ashore of ice water slurries through the static wrapping 
of the ice in a waterproof plastic barrier and the extraction 
of the fresh water melt. 

If the technology can be developed to transfer Antarctic 
ice at reasonable cost to the coasts of arid and semiarid 
regions, immense opportunities would accrue to governments 
with the resources to exploit the technology. A government 
like that of the United States could work in a cooperative, 
positive way with the government of the Republic of South 
Africa to open up approximately 50,000 Square miles of 
territory to agriculture and grazing in the more arid sections 
of the Cape of Good Hope and Orange Free State Provinces. 

The projects would be based ultimately on the technology 
necessary (1) to extract the maximum water from the existing 
rainfall and (2) to transfer and distribute additional water 
from Antarctic tabular bergs. The financial support for the 
water supply and agricultural project would involve invest- 
ment in a water distribution system and encouragement for 
settlement in the new territories. Such financing might be 
accomplished through a combination of U.S., Western European, 
South African government, World Bank, and private capital. 
The governments involved in cooperating in such a venture 
would be presented with opportunities to effect change in the 
hypothetical case of South Africa, without recourse to the 


hackneyed, aging and brutal techniques of insurgent and 
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Meeonter-insurgent warfare, without the advantages accruing 
in insurgency situations to the weaker and less developed 
Eurasian communist countries, and with the strong possibility 
of effective change and political accommodation based on 
rising standards of living within the Republic of South Africa. 
The United States, in particular, with its advanced technology, 
large financial resources, powerful and advanced allies, and 
traditional predominance at sea is faced with the opportunity 
Momuse its natural strengths to support successful political 
accommodation and settlement in Southern Africa. 

In other parts of the world, in situations less traumatic 
maa in South Africa, the U.S. has the opportunity through 
Bie exploitation of its technology to intensify agricultural 
activity in some areasS and open up agricultural and asso- 
Clated economic activity in others. On the west coast of 
North America, the California and U.S. Governments are faced 
memeiead three-year drought crisis (1975, 1976, 1977). The 
crisis may aernasee the positive supply of water in addition 
to the present negative provision of relief funds in order 
tO master the situation and continue to support the existing 
Overpopulation of California relative to scarce and ultimately 
limited water supplies. In the absence of the timely availa- 
bility of water from inter-region transfer and nuclear 
deslinization, Antarctic ice transfer probably will emerge 
as an imaginative, relatively low-cost alternative for the 


provision of fresh water to areas south of the equator. 


96 





With knowledge of the processes of melting, scientists, 
engineers, and governments could provide water to arid and 


semiarid regions north of the equator. 
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Nomex 


SOLUTION OF TRANSIENT THERMAL CONDUCTION EQUATION 


The comparison between the theoretical curves and 
experimental data points contained in Chapter VI, Figures 
100-105 were obtained utilizing the computer program, 
eeecuce Data on Temperatures Internal to Ice", Figure B-l. 
The terminology used for the text and computer programs is 
/meecea in Table B-I. 

The computer program, “Select a best S", Figure B-2, 
was utilized to assist in obtaining values for thermal 
diffusivity, ice regression rate, and Xe for program Figure 
Bt. The value ‘S" is found by subtracting the experimentally 
determined values from the calculated values, squaring the 
difference and summing over the set of experimental data 
Meants. The theoretical function, which provides calculated 
Values, was varied by changing the parameters a, Kor and ies 
A minimum value of "S" represented the best set of above 
values. 

After experimenting with the aforementioned parameters, 


a minimum value of "S" of 1329 was attained using ve = 2.2 Sor 


anda= Peo x10; these values of Ve and a are utilized 
miethe program shown in Figure B-l. 


The theory used to generate computer program 1 was the 


solution of the heat conduction equation as follows: 


ous 





Symbol in 
Computer Program 
me 9831 


A/A1L 


XO 


aN 


TQ 


X9 


W 


ieee 1 


NOMENCLATURE 
Symbol 
in 
Text Define ton 
H heat of fusion for ice 
q Detter tiux eo melt ice 
OL mmeiinal difrtusiyi ty 
of ice 
K EVe@mMcl  COnauUet dvl cy 
of ice 
Ve ice regression rate 
Ky distance from edge 
Gpeatece co first 
thermocouple 
Le surface temperature 
of ice 
T; interlor temperature 
of ice 
AX thermocouple spacing 
ii time intervals (1 - 6) 
0 density of ice 
W V_/o 


= 


Unies 


joule/kg 
p 
joule/m™ sec 


ao aXe 
jyoule/msec °K 
m/sec 

MmecEers 

7€ 

2 


meters 
minutes 
kaym- 


1l/meters 





An energy balance at the water-ice interface from above is: 





(ese) 


Of the heat transferred from the water the amount of heat 


which causes melting is: 


aX 
= aera = -? 
Q H 0 at HOV. sry) 


For the purpose of analysis, the temperatures within the 


ice are modeled using the one-dimensional unsteady heat 


transfer equation as follows: 
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The following assumptions were made: 
it) q, = constant 
2) Steady Peace Melting Occurs 
5) V, = constant = velocity of surface of ice 


The boundary conditions at the water-ice interface are as 


follows: 


oT 
ea) = a 
ox y=X = Hove qe (B-4) 
x 
T(X.,t] = T (B-5) 


[Mie DOUndary condition for the ice interior is: 


coe tele = De (B=6) 


which mathematically is as X + »~. 


Define a variable which involves both space and time 


Mete that the solution of eq. (B-3) will use € as an 
independent variable. Consequently the temperature profile 
propagates with the same shape into the ice as the ice 


melts. Combining equations, one can derive the following: 


a eeg the Bi = 


Low 





where the function "f" is given by: 
eet. = (T -T.] £12) (B-9) 


We now define: 


US oe (3-170) 
Vs ie 
Metabo r ' og Y 8 
V 
d Ss 7 
ae > O° 
V 
" 6} 
-V,, §/0 
yp = ae 12 18 


Taking the derivative of f£f with respect to &€ and then 


integrating to solve for £f as follows: 


ae 
arp = ae 
SVE / c Nee erect 
a oe ae c=) a eae +b 
QL 
c. =saiame eae: 
av. 570 
P= ys Ce + b (B=11) 
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Mepitzing boundary condition (6) and combining eq. (B-9) 


Sma (B-1l) we obtain the following: 


eae. > (T.-T.] [ce © +b] 


Memttzing boundary condition (5) when X 


following: 


It therefore follows that: 


TX, ,t] = aie - (ya 


ed = ieee eey | hoY 7c) = VoE/o 
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= T. (e172) 


X we obtain the 


0 

(B-13) 
eine) fear Vel 
[L - ae 

(B-14) 


Ls 


d 
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ms 


I 
! 





We note that when A Xa € = 0, therefore: 





ceeeeetiow es K7a V_ [T_-T. | 


The final result is as follows: 


(We 
Sho a V , [Ho eh ae (B-16) 
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Regimes of flow and wave producing conditions as 
functions of Reynolds and Proude numbers. 


Regimes of flow and wave producing conditions as 
functions of Reynolds and Froude Numbers. 
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2. NPS water tunnel, 
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8. Effectiveness of heat transfer. 
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19. Scale used to weigh ice and container. 
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GLACIAL ICE THERMOCOUPLE POSITIONS 
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4la. Thermocouple arrangement; Test 4. 
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Thermocouple with digital readout. 
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GLACIAL ICE THERMOCOUPLE POSITIONS 
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Forward section separation; Test 6. 
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